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Looking at the Gregory-Laflamme instability through quasi-normal modes 
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We study evolution of gravitational perturbations of black strings. It is well known that for 
all wavenumber less than some threshold value, the black string is unstable against scalar type 
of gravitational perturbations, which is named the Gregory-Laflamme instability. Using numerical 
methods, we find the quasinormal modes and time-domain profiles of the black string perturbations 
in the stable sector and also show the appearance of the Gregory-Laflamme instability in the time 
domain. The dependence of the black string quasinormal spectrum and late time tails on such 
parameters as the wave vector and the number of extra dimensions is discussed. There is a numer- 
ical evidence that in the threshold point of instability the static solution of the wave equation is 
dominant. For wavenumbers slightly larger than the threshold value, in the region of stability, we 
see tiny oscillations with very small damping rate. While, for wavenumbers slightly smaller than 
the threshold value, in the region of the Gregory-Laflamme instability, we observe tiny oscillations 
with very small growth rate. We also find the level crossing of imaginary part of quasinormal modes 
between the fundamental mode and the first overtone mode, which accounts for the peculiar time 
domain profiles. 
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I. INTRODUCTION 



Unlike four dimensional Einstein gravity, which allows 
existence of black holes, higher dimensional theories, such 
as the brane-world scenarios and string theory, allow ex- 
istence of a number of "black" objects: higher dimen- 
sional black holes, black strings and branes, black rings 
and saturns and others. In higher than four dimensions 
we lack the uniqueness theorem, so that stability may 
be the criteria which will select physical solutions among 
this variety of solutions. Up to now, we know that higher 
dimensional Reissner-Nordstrom-de Sitter black holes are 
stable pi in the Einstein gravity. On the contrary, black 
holes in Gauss-Bonnet (GB) gravity are unstable for large 
GB coupling for D = 5, 6 |2|, where D is the total number 
of space-time dimensions. Kaluza-Klein black holes with 
squashed horizon are stable against lowest zero mode 
perturbations [3|. Unlike, Kaluza-Klein black holes, the 
black string metric is a solution of the Einstein equations 
in five or higher dimensional gravity that has a factor- 
ized form consisting of the Tangherlini black hole and an 
extra flat dimension [6J. According to the brane-world 
scenarios, if the matter localised on the brane undergoes 
gravitational collapse, a black hole with the horizon ex- 
tended to the transverse extra direction will form. This 
object looks like a black hole on the brane, but is, in fact, 
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a black string in the full D-dimensional theory. 

It is well known that such black strings suffer from the 
so-called Gregory-Laflamme instability, which is the long 
wavelength gravitational instability of the scalar type of 
the metric perturbations (4j, [fj]. The Gregory-Laflamme 
instability has been intensively studied during the recent 
decade [fj| and the threshold values of the wave vector k 
at which the instability appears are known [7|. In the 
present paper we are aimed at studying the evolution 
of linear perturbations of D-dimensional black strings 
in time and frequency domains. This task is motivated 
mainly by the two reasons: first to realize what happens 
on the edge of instability of black strings and how the 
perturbations will develop in time. Second, the so-called 
quasinormal modes of a stable black string might be an 
observational characteristic for the future Large Hadron 
Collider Experiments, if such objects as black strings ex- 
ist. 

The latter needs some more explanation. In this re- 
search, we shall show that, if a stable D-dimensional 
black string is gravitationally perturbed, it will undergo 
damped oscillations, called quasinormal ringing, similar 
to that of a black hole Q. At asymptotically late time, 
this quasinormal ringing goes over into the power-law 
tails. The quasinormal modes and asymptotic tails are 
very well studied for D-dimensional black holes [lj| and 
for black holes localised on the brane [11]. We find here 
that the quasinormal ringing of black strings has a num- 
ber of differences from that of D-dimensional black holes 
[lfj |. especially when approaching the edge of instabil- 
ity. In particular, we find a numerical evidence that in 
the threshold point of instability, the static solution of 
the wave equation is dominant. In the region of the sta- 



bility, for k slightly larger than the threshold value, we 
see modes with tiny oscillation frequencies and damping 
rates. In the region of the Gregory-Laflamme instability 
for k slightly smaller than threshold value, we observe 
tiny oscillations and instability growth. We also find the 
level crossing of imaginary part of quasinormal modes 
between the fundamental mode and the first overtone 
mode. 

The paper is organized as follows. Sec II gives ba- 
sic formulas for the black string perturbations and the 
wave equation for the scalar type of gravitational per- 
turbations. Sec III analyses quasinormal modes of black 
strings in frequency domain with the help of the Frobe- 
nius method, and in time domain by the Gundlah-Price- 
Pullin method. We also discuss features of the Gregory- 
Laflamme instability using the obtained results. The fi- 
nal section is devoted to the conclusion. 



II. BASIC FORMULAE 

In this section, we shall briefly review the results of 
the paper [3], where the wave equation for the scalar 
type of gravitational perturbations were obtained. This 
wave equation will be our starting point for numerical 
investigation. 

For the static string in D = n + 4 space-time dimen- 
sions, the background metric can be written as 

dr 2 



ds z = g^dx^dx" = -f{r)dt z + — + r 2 dQ, 2 n+1 + dz l 



where 



(1) 



fir) 



1 



and df^ +1 is the metric on a unit (n+ l)-sphere. Various 
properties of black string have been extensively studied 
in recent years and we refer a reader to the papers [1 21 ] — 
[27| for more detailed information on black strings. 

The z-direction is periodically identified by the relation 
z = z + 2nR. We study perturbations of an (n + 1)- 
spherically symmetric solution with the Killing vector in 
z-direction. Therefore, we can write perturbations in the 
form 



5g^ 



r j.kz 



a^v(t,r), 



R 



771 G Z. 



The perturbed vacuum Einstein equations have the form 



SRuv = 



(2) 



The perturbations can be reduced to the form, where the 
only non-vanishing components of a^ v are 



a-tt = h t 



(Xtr 



/i r , Oj ZZ n z , 



Ifcfly . 



The linearized Einstein equations give a set of cou- 
pled equations determining the four radial profiles above. 
However, we may eliminate h v , h r and h t from these 
equations in order to produce a single second order equa- 
tion for h z : 



h z = f{r) 2 K+p(r)ti z +q{r)h z 



(3) 



where 



p(r) 



f{rf 



1 



4(2 + n)k 2 r 2 



q(r) = -k z f{r) 
Defining 



\ /(r) 2k 2 r 2 + n(n + l)(r + /r)« 

2k 2 r 2 -n(n + 3)(r + /r) n 



2k 2 r 2 +n(n + l)(r+/r) n 



h z (t,r) 



r -(n-l)/2 

2k 2 r 2 +n(n + l)(r+/r) n v ' h 



we can reduce the equation ^ to the wave-like equation 

d 2 d 2 



dt 2 dr 2 

dr 



+ V{r) * = 



(4) 



where dr+ — -rrp- is the tortoise coordinate. Here, the 
effective potential V(r) is given by 

V(r) = ™ ^ 

4r 2 (2k 2 r 2 +n{n+l)(r+/r) n ) 2 ' 

where 

U(r) = 16k 6 r 6 + 4:k\ 4 (n + 5)(3/(r) - 2n + 3n/(r)) - 
-4k 2 r 2 n(n + l) (n(n+5) + /(r)(2n 2 + 7n+9)) 

r) ' 
The above effective potential does not vanish at asymp- 
totic infinity but has an effective " mass" term, containing 
k, at the spatial infinity. 



III. EVOLUTION OF PERTURBATIONS 

ANALYZED WITH THE FROBENIUS METHOD 

AND TIME DOMAIN INTEGRATION 

TECHNIQUE 

First of all, let us briefly describe the two methods 
which we used here: the Frobenius method (frequency 
domain) and the Gundlach-Price-Pullin method (time 
domain) . 

In time domain, we study the black string ringing using 
a numerical characteristic integration method [28], that 
uses the light-cone variables u — i — r* and v = t + r* . In 
the characteristic initial value problem, initial data are 
specified on the two null surfaces u = uq and v = vq. 
The discretization scheme we used, is 



*(7V) = ®(W) + V(E) - *(S) 



(5) 



2 V{W)^(W) + V{E)^(E) 



0(A 4 ) 



TABLE I: Fundamental mode (wo) found by the time-domain 
integration, first (ui) and second (0^2) overtones of spherically 
symmetric black string perturbations (n = 1) found by the 
Frobenius method. As k grows, the first overtone decreases 
its imaginary part, becoming the fundamental mode, reaching 
purely real frequency (quasi-resonance) at k » 0.94 and then 
disappearing. The second overtone, as well as the fundamen- 
tal purely imaginary mode, increases its damping rate with 
k. 



k 


io (t-d) 


wi (Frob.) 


0.84 


+0.01H 


0.802 - 0.0173; 


0.85 


+0.008* 


0.810 - 0.0157; 


0.86 


+0.005; 


0.818 - 0.0140; 


0.87 


+0.002; 


0.827 - 0.0124; 


0.88 


-O.OOli 


0.835 - 0.0107; 


0.89 


-0.004; 


0.844 - 0.0090; 


0.90 


-0.007; 


0.852 - 0.0074; 


0.91 


-0.011? 


0.861 - 0.0058i 


0.92 


-0.014i 


0.870 - 0.0042i 


0.93 


-0.0174 


0.879 - 0.0025; 


0.94 


-0.021* 


0.888 - 0.0007; 



k 


u>2 (Frob.) 


0.8 


0.361 - 0.632; 


0.9 


0.373 - 0.679; 


1.0 


0.378 - 0.724* 


1.1 


0.374 - 0.767; 


1.2 


0.364 - 0.806; 


1.3 


0.347 - 0.841; 


1.4 


0.324 - 0.870; 


1.5 


0.294 - 0.894; 



where we have used the following definitions for the 
points: N = (u+A, v+A), W = {u+A, v), E = (u, v+A) 
and S — (u, v). 

In frequency domain we used the well-known Frobe- 
nius method [29(. In order to study the QN spectrum 
in frequency domain, we have separated time and radial 
coordinates in ([3]) 



h z (t,r) 



'Mr). 



Here h u (r) satisfies the quasinormal mode boundary con- 
ditions, which are purely ingoing wave at the event hori- 
zon and purely outgoing wave at the spatial infinity. 
Thus, the appropriate Frobenius series are 



i(2u> 



y. \ —vjJT+jn 



where a = 



r 
1 k 2 ) 



e iVu> 2 -k 2 r r (n+3+a)/2 / \ 



(6) 



Vu' 



for n = \ and a = for n > 1. 



j 2 -k 2 

It is crucial that y(r) must be regular at the event horizon 
and at the spatial infinity and can be expanded as 



*r)=£>(i-^y 



i=0 

Substituting © into §S§, we have found that the coeffi- 
cients cii satisfy a (3n + 5)-term recurrence relation. We 
found the coefficients of the recurrence relation, and then 
we obtained the equation with the infinite continued frac- 
tion, which is algebraic equation with respect to the QN 
frequency u>. Numerical solutions of this algebraic equa- 
tion give us the QN spectrum. 

Now let us discuss the obtained results for QN modes 
and time domain profiles. The Frobenius method for the 
considered cumbersome potential gives rise to a num- 
ber of technical difficulties: first, the convergence of the 



TABLE II: Dominant QNMs of spherically symmetric black 
string perturbations (n = 2) found with the time-domain and 
the Frobenius method. For k > 1.4 the wo mode is not domi- 
nating anymore and is difficult for detection by the time do- 
main integration. 



k 


wo (t-d) 


wi(t-d) 


u>i (Frobenius) 


1.3 


-0.012; 


1.418-0.14H 


1.396-0.118; 


1.4 


-0.047; 


1.491 - 0.099J 


1.483-0.107; 


1.5 


n/a 


1.584 - 0.095J 


1.570 - 0.096; 


1.6 


n/a 


1.674 - 0.088J 


1.659 - 0.086; 


1.7 


n/a 


1.763 - 0.081; 


n/a 


1.8 


n/a 


1.852 - 0.077; 


n/a 


1.9 


n/a 


1.936 - 0.054; 


n/a 


2.0 


n/a 


2.032 - 0.046J 


n/a 



TABLE III: Dominant QNMs of spherically symmetric black 
string perturbations (n — 3) found with the time-domain 
method. For k > 2 the u>i mode is a dominating one. 



k 


uoa 


LOl 


1.6 


-0.007; 


1.869 - 0.214; 


1.7 


-0.043; 


1.930-0.211; 


1.8 


-0.082; 


2.013-0.212; 


1.9 


-0.130; 


2.106-0.185; 


2.0 


-0.250; 


2.223 - 0.179; 


2.1 


n/a 


2.309-0.172; 


2.2 


n/a 


2.400 - 0.167; 


2.3 


n/a 


2.472 -0.161; 



Frobenius series is rather slow. Second, when searching 
for the solutions of the algebraic equation in the region 
close to the threshold of instability, one needs very good 
initial guess for u> to "fall" into the minimum of the con- 
tinued fraction equation. This can be easily understood. 
As we shall see from the time domain integration, the 
dominant solution in the threshold point corresponds to 
some static solution with u> = 0, so that nearby funda- 
mental mode has tiny real and imaginary parts. The 
Frobenius method naturally requires slow convergence 
and excellent initial guess for to for such small u>. 

In tables I, we have listed fundamental mode ujq, first 
lo\ and second tt>2 overtone of spherically symmetric black 
string perturbations in the case of n = 1. First, we notice 
the level crossing of imaginary part of quasinormal modes 
between the fundamental mode and the first overtone 
mode. This level crossing is peculiar to the black strings. 
As k grows, the first overtone u\ decreases its imaginary 
part, becoming the fundamental mode, reaching purely 
real frequency (quasi-resonance) at k « 0.94 and then 
disappearing. The second overtone, as well as the fun- 
damental purely imaginary mode, increases its damping 
rate with k. 

In table II, we listed dominant QNMs of spherically 
symmetric black string perturbations in the case of n = 2. 
We can see that the fundamental mode near the threshold 
of instability has no real part. Yet, higher overtones have 
detectable real part. 

In table III, we listed dominant QNMs of spherically 
symmetric black string perturbations in the case of n = 



Re(w) 




Iirfw) 



-0.8 
-1.0 



FIG. 1: Real and imaginary part of first two overtones for n = 2 black string perturbations as a function of k. The first overtone 
wi (blue) approaches lu = k (red line). The second overtone u>2 (green) becomes pure imaginary in stable region (k ~ 1). 
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FIG. 2: Time-domain profiles of black string perturbations 
for k = 2.5 n = 2 (red, top), n — 3 (orange), n = 4 (green), 
n — 5 (blue, bottom). Late-time decay of perturbations for 
n > 3 is oc t~(' I+6 )/ 2 . For lower n the law of decay is different: 
oc t -0 ' 93 for n = 1 and oc i" 1 ' 2 for n = 2. 



3. We can see the level crossing of imaginary part of 
quasinormal modes between the fundamental mode and 
the first overtone mode. 

The real part of the second mode (wi in tables I, III) 
asymptotes to k at large k, while the imaginary part 
monotonically decreases when k is increasing. The third 
mode does not asymptotes to k, but has monotonically 
decreasing real and imaginary parts as can be seen from 
Fig. 1. Because of the risk to "fall" into another over- 
tone, in order to obtain the higher overtones, we had to 
start from the D-dimensional black holes with k = 0, 
for which the QN frequencies are known , and then to 
"move" towards higher k in the Frobenius method with 
a very small step (see Fig. 1). 

An essential advantage of the time domain method in 
comparison with the Frobenius method is that we do not 
have any decreasing of the convergence or loss of accuracy 
when approaching the point of instability. Therefore our 



time domain method is more complete than the frequency 
one, at least for the dominant mode, which can always 
be extracted from the time domain picture. In Fig. 2, 
one can see examples of time domain profiles for various 
n and a fixed k. There one can see that for n > 3 the 
intermediate late time asymptotic is power law like 

# oc r (n+6)/2 , n>3, 

while for other n, the asymptotics are 

•Foci" - 93 , for n = l, 



* oc f 



for n = 2. 



Let us note that this asymptotic apparently should be 
considered as intermediate. They are expected to go over 
into other power law ones at very late times, as it takes 
place for massive fields in general [301, teMi uSi [2JI- 

Let us note, that k plays the role of the effective 
mass. At asymptotically late time we observe power- 
law damped tails, which have oscillation frequency equal 
to k, resembling asymptotical behavior of massive fields 
near Schwarzschild black holes. The first overtone's 
behavior is qualitatively similar to that of the funda- 
mental mode for massive fields of higher-dimensional 
Schwarzschild black holes [32|]: we can see long-lived os- 
cillations, which can be infinitely long lived modes called 
the quasi- resonances [33|. The analytical explanation of 
existence of the quasi-resonances was found in [34J . 

• For D = 5 (n = 1), as k grows, the imaginary part 
of the first overtone quickly decreases and vanishes 
for some threshold value of k, while its real part 
stays smaller than the threshold value (see Table 
|T]). After the threshold value of k is reached, the 
first overtone "disappears". 

• For D > 6 (n > 2), the imaginary part of the first 
overtone becomes small for large k, while the real 
part asymptotically approaches k (see Fig. [1]). 



Even though the first overtone of the spherically sym- 
metric black strings behaves similarly to the funda- 
mental mode of massive fields near higher-dimensional 
Schwarzschild black holes, the other modes have com- 
pletely different behavior. The fundamental mode of a 
black string perturbation is purely imaginary. It grows 
for small values of k, leading to instability of the black 
string. This behavior is common for the unstable modes. 

Indeed, let us multiply the equation (|4|) by the complex 
conjugated function \&* and assume that the dependance 



on time is ^(t, r*) 



\I/(?v). Let us study the inte- 



gral of the obtained equation 



**(r*) 



dr 2 . 



+ uj 2 \y(r,)\ 2 -V\y(r+)\ 2 )dr+ 



Integration of the first term by parts gives 



I = tf*(r*) 




dr+ 



Lo 2 \V(n)\ 2 ~V\t>(n)\ 2 



d^{r+) 



dr* 



dr* = 0. 



Taking into account the boundary conditions © , we find 
that the imaginary part of the integral is 



Im(I) = Re{y/oj 2 - fc 2 )|*(oo)| 2 + i?e(w)|*(-oo)| 2 + 

OO 

+2Re(uj)Im{uj) f |*(r*)| 2 dr* = 0. 



Since the sign of Re(Vu> 2 — k 2 ) coincides with the sign 
of i?e(w), the non-zero real part of the quasi- normal 
frequency implies that the imaginary part is negative. 
Therefore, the unstable modes (Im(ui) > 0) must have 
zero real part. In other words, unstable modes cannot be 
oscillating. 

In the stability region, the fundamental mode is also 
purely imaginary, but damped mode, whose damping 
rate grows quickly with k. Because of its quick growth, 
this mode cannot be considered as the fundamental one 
for large k. In fact, for larger k (see Tables I-III), the 
first overtone turns out to be the fundamental mode (the 
mode with the largest lifetime). The real part of the 
second overtone decreases, as k grows, and reaches zero 
for some value of k, while the imaginary part remains 
negative. 

Let us now look at the Fig. 3. At moderately large 
values of k, sufficiently far from instability, the profile 
has the same form as that for massive fields, yet, when 
approaching the instability point the real oscillation fre- 
quency (Re(u>)) and the decay rate (Im(u>) < 0) decrease 
considerably. After crossing the instability point we ob- 
serve, that starting from some tiny values, Im(uj) > are 
slowly increasing (while Re(uj) is still zero for the funda- 
mental mode and tiny oscillations, observed in the time 
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FIG. 3: Time-domain profiles of black string perturbations 
for ri — 1 k — 0.84 (magenta, top), k = 0.87 (red), k = 
0.88 (orange), k — 0.9 (green), k = 1.1 (blue, bottom). We 
can see two concurrent modes: for large k the oscillating one 
dominates , near the critical value of k the dominant mode 
does not oscillate (looks like exponential tail), for unstable 
values of k the dominant mode grows. The plot is logarithmic, 
so that straight lines correspond to an exponential decay. 



domain, come from the next decayed mode). Therefore 
we conclude, that the there is some static solution to = 
of the wave equation (4), which shows itself exactly in 
the threshold point of instability. We would say that this 
picture of instability is natural, if the instability develops 
on the fundamental mode. However, if instability occurs 
at higher multipoles £, as it takes place for instance in 
the Gauss-Bonnet theory [2| , the picture of instability is 
quite different: growing modes appear only after rather 
long period of decaying oscillations (see [2|]). Note also, 
that here we confirmed the threshold values of k found 
in [7| with a very good accuracy by the time domain in- 
tegration (see for instance Fig. 3 for n = 1). Thus the 
threshold values are: k — 0.876 for n — 1, k = 1.269 for 
n = 2, k = 1.581 for n = 3, k = 1.849 for n = 4. 

Finally, in Fig. 4, we can see the region of the profiles 
where the period of the quasinormal ringing goes over 
into the power-law tail behavior. Close to the critical 
point, there exists a period where the oscillation ceases. 
This is because the pure damping mode becomes the fun- 
damental mode near the critical point. 



IV. CONCLUSION 

We have numerically studied the Gregory-Laflammc 
instability through quasinormal modes. Let us stress 
three main results obtained here: 

1) We have found the quasinormal modes and late- 
time tails for scalar type of gravitational perturbations 
of D-dimensional black strings for various D, that is for 
the type of perturbations where the Gregory-Laflammc 
instability forms in the long wavelength regime. 

2) The time domain profiles indicate that the threshold 




FIG. 4: Time-domain profiles of black string perturbations for n — 4 k — 2.2 (red, left top), k = 2.3 (orange, right top), k = 2.4 
(green, left bottom), k = 2.5 (blue, right bottom). Quasi-normal ringing and tails have the same frequency of oscillation which 
is close to k. One can see a period where the oscillation ceases close to the critical point. 



instability value of k corresponds to dominance of some 
static solution lo = 0. 

3) Near the instability point (in k) the fundamental 
mode is pure imaginary (non-oscillating), and, as k is in- 
creasing, the lifetime of the second mode is increasing, 
so that at some moderate k both modes are dominating 
at the late time of the ringing. At larger k, the domi- 
nance goes over to the second (oscillating) mode, as to 
the longer lived one. 

Our research could be improved in a number of ways. 
First of all, one could compute QNMs for higher multi- 
pole numbers, starting from the effective potential de- 
rived in [35| and also for other types of gravitational 
perturbation. Though vector and tensor types of per- 
turbations do not contain instabilities, such investigation 
would give us complete data on QNMs and evolution of 
gravitational perturbations. 

The main limitation of our analysis is that we cannot 



say what happen with unstable black string some time 
since the moment of initial perturbations: the pertur- 
bations will grow and become large, so that the linear 
approximation will not be valid anymore. However, it is 
beyond the scope of this paper. 

Acknowledgments 



R. A. K. was supported by the Japan Society for the 
Promotion of Science (JSPS), Japan. K.M. was sup- 
ported by JSPS Grant-in- Aid for Scientific Research No. 
19 • 3715. J.S. was supported by the Japan-U.K. Re- 
search Cooperative Program, Grant-in- Aid for Scientific 
Research Fund of the Ministry of Education, Science and 
Culture of Japan No. 18540262 and No. 17340075. A. Z. 
was supported by Fundagao de Amparo a Pesquisa do 
Estado de Sao Paulo (FAPESP), Brazil. 



[1] H. Kodama and A. Ishibashi, Prog. Theor. Phys. 
Ill, 29 (2004) |arXiv:he p-th/0308128 ; A. Ishibashi 
and H. Kodama, Prog. Theor. Phys. 110, 901 



(2003) [arXiv:hep-th/0305185| ; R. A. Konoplya and 
A. Zhidenko, Nucl. Phys. B 777, 182 (2007) 



arXiv:hep-th/0703231 



[2] R. A. Konoplya and A. Zhidenko, Phys. Rev. D 77, 

104004 (2008) larXiv :0802.0267 [hep-th]; M. Beroiz, 
G. Dot ti and R. J. Gleiser, P hys. Rev. D 76, 024012 

(2007) [arXiv:hep- th/0703074 . 

[3] H. Ishihara, M. Kimura, R. A. Konoplya, K. Murata, 
J. Sod a and A. Zhidenko, Phys. Rev. D 77, 084019 

(2008) arXiv:0802.0655 [hep-th]]; M. Kimura, K. Mu- 
rata, H. Ishihara and J. Soda, Phys. Rev. D 77, 064015 
(2008) arXiv:0712.4202 [hep-th]]. 

[4] R. Gre gory and R. Laf lamme, Phys. Rev. Lett. 70, 2837 

(1993) [arXiv:hep- th/9301052] . 

[5] R. Gre gory and R. La flamme, Nucl. Phys. B 428, 399 

(1994) [arXiv:hep-th/94040711 . 

[6] T. Harmark, V. Niarchos and N. A. Obers, Class. Quant. 
Grav. 24, Rl (2007) [arXiv:hep-th/0701022] . 

[7] J. L. H ovdebo and R. C. Myers , Phys. Rev. D 73, 084013 
(2006) arXiv:hep-th/0601079 . 

[8] K. Kokkotas, B. Schmidt, Liv. Rev. Rel. 02, 2, 1999 

[9] R. A. Konoplya, Ph ys. Rev. D 68, 124017 (2003) 
arXiv:hep-th/0309030 . 
[10] R. Konoplya, Phys Rev. D 71, 024038 (2005) 
[arXiv:hep-th/0410057| ; V. Cardoso, J. P. S. Lemos 
and S. Yoshida, JHEP 0312 (2003) 041 
|arXiv:hep-t h/0311260 ; H. T. Cho, A. S. Cornell, 
J. Dou kas and W. Naylor, P hys. Rev. D 75 (2007) 

104005 arXiv:hep-th/0701193 ; R. G. Daghigh, G. Kun- 
statter and J. Ziprick, Class. Quant. Grav. 24 (2007) 
1981 [arXiv:gr-qc/0611139| ; F. W. Shu and Y G. Shen , 
JHEP 0608 (2006) 087 [arXiv:hep-th /0605128 ; 
R. A. Konoplya, P hys. Rev. D 68, 024018 (2003) 

arXiv:gr-qc/0303052 ; R. A. Konoplya, Phys. Rev. D 66, 
044009 (2002) arXiv:hep-th/0205142 ; R. G. Daghigh, 
G. Kunstatter, D. Ostapchuk and V. Bagnulo, Class. 
Quant. Grav. 23 (2006) 5101 [arXiv:gr-qc /0604073 ; 
R. A. Konoplya and E. Abdalla , Phy s. Rev. D 71, 
084015 (2005) [arXiv:hep-th/ 0503029 ; E. Abdalla, 
R. A. Konopl ya and C. Molina, P hys. Rev. D 72, 
084006 (2005) arXiv:hep-th/0507100 ; A. Lopez-Ortega, 
arXiv:0706.2933 [gr-qcj; C. Moliha T Phys. Rev. D 68 
(2003) 064007 arXiv:gr-qc/0304053 . 

[11] P. Kant i and R. A. K onoplya, Phys. Rev. D 73, 044002 
(2006) arXi v:hep-th/0512257j ; P. Kanti, R. A. Kono- 
plya and A. Zhid enko, Phys. Rev. D 74, 064008 (2006) 
arXiv:gr-qc/0607048 ; M. Nozawa and T. Kobayashi, 
larXiv:0803.3317l [hep-th]; U. A. al-Binn i and G. Siopsis, 
Phys. Rev. D 76 (2007) 104031 [arXiv:0708.3"363l [hep- 
th]]; , 

[12] L. h. Liu and B. Wang, arXiv:0803.0455 [hep-th]. 



[13 
[14 

[is; 
[is; 

[17; 
[is; 

[19 

[20 
[21 

[22 

[23 

[24 

[25 

[26 

[27 
[28 

[29 
[30 

[31 

[32 
[33 
[34 



Y. Kurita and H. Ishihara, Class. Quant. Grav. 25, 
085006 (2008) arXiv:0801.2842 [hep-th]]. 
Y. Brihaye, T. Delsate and E . Radu, Phys. Lett. B 662, 
264 (2008) [arXiv:0710.4034l [hep-th]]. [lj| 

B. Kleihaus and J. Kunz. larXTv:0710.1726l [hep-th], [j]| 
A. Bernamonti, M. M. Caldarelli, D. Klemm, R. Olea, 

C. Sieg an d E. Zorzan, JHEP 0801, 061 (2008) 
arXiv:0708.2402 [hep-th]]. 

Y. Brihaye and T. Delsate, Class. Quant. Grav. 24, 4691 

(2007) [arXiv:gr-qc/0 703146 . 

K. Murata, J. Soda and S. Kanno, Phys. Rev. D 75, 

104017 (2007) arXiv:gr-qc/0703029; ; K. Murata, J. Soda 

and S. Kanno, arXiv:gr-qc/0701137 

Y. Brihaye, B. Hartmann and E. Radu, Phys. Rev. D 72, 

104008 (2005) [arXiv:hep-th/0508028l . 

A. Tomimatsu, Phy s. Rev. D 71, 124044 (2005) 
[arXiv:gr- qc/0506120\ 

T. Kobayashi and T. Tanaka, Phys. Rev. D 71, 084005 
(2005) [arXiv:gr-qc70 412139 . 

V. Car doso and J. P. S. Lemo s, Phys. Lett. B 621, 219 
(2005) [arXiv:hep-th/0412078] . [H 

B. Kol, Phys. Rept. 422, 119 (2006) 
[arXiv:hep-th70 411240 . 

M. I. Park, C lass. Quant. Grav. 22, 2607 (2005) 

[arXiv:hep- th/0405045 . 

T. Tamaki, S. Kanno and J. Soda, Phys. Rev. D 69, 



024010 (2004) arXiv:hep-th/0307278 . 

S. Kanno and J. Soda, Class. Quant. Grav. 21, 



1915 



(2004) [arXiv:gr-qc/ 03 11074 . 



S. Kanno and J. Soda, arXiv:gr-qc/0312107 
C. Gundlach, R. H. Price and J. Pullin, Phys. Rev. D 49, 
883 (1994) [arXiv:gr-qc/9307009j . 

E. W. Leaver, Proc. Roy. Soc. Lond. A 402, 285 (1985). 
L. E. Simone and C. M. Will, Class. Quant. Grav. 9 
(1992) 963. 

R. A. Konoplya, P hys. Lett. B 550, 117 (2002) 
arXiv:gr-qc/0210105 

A. Zhidenko, Phys . Rev. D 74, 064017 (2006) 
arXiv:gr-qc/0607133 . 

A. Ohashi and M. a. Sakagami, Class. Quant. Grav. 21 
(2004) 3973 [arXiv:g r-qc/0407009 

R. A. Konoplya and A. V~~Z hidenko, Phys. Lett. B 
609, 377 (2005) a rXiv:gr-qc/041 1059 ; R. A. Kono- 
plya and A. Zhidenko, Phys. Rev. D 73, 124040 (2006) 
[arXiv:gr-qc/06050l3] . 
[35] H. Kudoh, Phys. Rev. D 73, 104034 (2006) 
[arXiv:hep-th/060200T] . 



